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We present results from 43 GHz (VLBA, six epochs from 2012.2 to 2013.2) and 86 GHz (GMVA, 
one epoch in 2012.4) observations toward the basis of the jet in the TeV Blazar Mrk 501. The 
43-GHz data analysis reveals a new feature located northeast of the radio core, with a flux density 
of several tens of mJy, perpendicularly to the jet axis. The 86-GHz image shows the jet feature 
located 0.75 mas southeast of the radio core, which is consistent with the previous result. The 
location of Gaussian model for 0.75 mas feature does not coincide with those for the jets in the 
43-GHz image, however, a distribution of emission is found. We also discuss the spectral indices 
of the core, the northeast feature, and the jet feature between 43 GHz and 86 GHz, which show 
flat-to-steep, steep, and flat-to-invert, respectively. 
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1. Introduction 

The innermost structure of relativistic jets in active galactic nuclei are considered to be related 
to the properties of the regions where the jets are formed, collimated, and accelerated. Ultra-high 
resolution Very Long Baseline Interferometry (VLBI) observations are the best tools to explore the 
parsec scale structure of jets. Recently, the Event Horizon Telescope at 230 GHz (1.3 mm) has 
revealed the off-axis jet structure within 0.3 mas from the core in two quasars [1, 2]. However, in 
BL Lac Objects, such off-axis jet structure have not been studied. 

The TeV blazar Mrk 501 is one of the best BL Lac Objects to resolve the innermost jet because 
of its proximity (z = 0.034, 1 mas = 0.66 pc) and brightness. The VLBI images show different jet 
position angles at different scales [3]. The previous Very Long Baseline Array (VLBA) 43 GHz and 
Global mm-VLBI Array (GMVA) 86 GHz images show clear limb-bright structure [4, 5], while 
they do not show clear off-axis jet structure at the innermost part. 

2. Observations and Data Reduction 

We observed this target in six different epochs by the VLBA at 43 GHz and one epoch by the 
GMVA at 86 GHz, with 8 IFs of 16 MHz bandwidth each. The total on-source time for Mrk 501 
was approximately one hour for the 43 GHz data, and four hours for the 86 GHz data. Phase 
and delay offset between different sub-bands were solved by using the calibrator 3C 345. For the 
43 GHz data, after the standard initial calibrations in AIPS, the final images were produced after 
iterations of CFEAN, phase and amplitude self-calibration processes in difmap. For the 86 GHz 
data, fringe fitting on Mrk 501 were performed with a small search window of a few tens of nsec 
and mHz in delay and rate respectively, and a low SNR threshold, after applying the delay and rate 
solutions of 3C 345. Fringe solutions of Mrk 501 were not derived among transatlantic baselines, 
however, weak coherent phases were found. Gaussian model fitting and phase self-calibration were 
done in difmap, without amplitude self-calibration due to low SNR of the visibilities. 

3. Results 

In Figure [l], we find a new component (labeled NE) located almost perpendicularly to the 
main jet axis (labeled Cl, C2, C3). The location (angular separation, position angle, and flux) 
of NE changes randomly. We also confirmed the existence of NE using the slice profiles of the 
CLEAN images across the center of Gaussian models for core and NE. There is a clear emission 
at the northeast of the core for all epochs with the peak flux of 36 mJy beam -1 , the position of 
0.18 mas from the core, and the FWHM of 0.24 mas in average. These results suggest that the 
limb-brightened structure of Mrk 501 apparently stars in the nearest vicinity of the core due to the 
small viewing angle. 

Figure 2 -Left shows the total-intensity image obtained by the GMVA at 86 GHz at a resolution 
of 205 jUas x 37 /ras (PA —5.11°) with fitted Gaussian models. We find a tentative jet knot located 
at (r, 6) = (0.75 mas, 159°) with a SNR of ~ 5(7. This knot could correspond to one detected 
by the GMVA at 86 GHz in [4], located at (r, 6) = (0.73 mas, 172°). The flux of the core is 
197 ± 59 mJy and that of the tentative jet knot is 19 ± 8 mJy. 
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Figure 1: 43 GHz uniform weighted VLBA images with fitted circular Gaussian components. The restoring 
beam is 0.39 mas x 0.14 mas in PA — 14°, plotted in the bottom-left corner. Observing date and the first 
contour (f.c.) are plotted above each map. The first contour are all set to three times the rms noise level of 
each map, increasing by a factor of 2. 

Since the 86-GHz data and the third-epoch VLBA data at 43 GHz were obtained within two 
weeks, we compare these two images. In the 43-GHz image, there are no Gaussian models at 
~ 0.75 mas southeast of the core, however, a distribution of emission is found. The spectral index 
of the integrated core fluxes over Gaussian models between 43 GHz and 86 GHz is a = —0.8 ± 0.5, 
where we define the flux density at observing frequency v as 5 V « v“. This flat-to-steep spectral 
tendency is consistent with the power law index of the core reported in [4] (a ~ —0.5 above 8 GHz). 
Using the model-fit flux of the 43-GHz image and the 3a upper limit in the 86-GHz image, the 
spectral index of NE is estimated to be a < —0.8, which is steeper than that for the radio core. We 
also estimate the spectral index at the southeast jet as a = —0.1 ±0.5 by integrating the flux within 
the same region between the two images using the AIPS task IMSTAT. Due to the large spectral 
index error, this is not a stringent result. However, despite of the steep spectrum of the core, the 
southeast jet spectrum is flat-to-inverted and it is inconsistent with one at lower frequencies [e.g., 
6]. This tendency is also seen in the spectral index map shown in Figure 2 -Right. 

4. Summary and future prospects 

We have explored the parsec-scale regions in the jet of Mrk 501 using the VLBA at 43 GHz 
and the GMVA at 86 GHz. We newly find the off-axis jet structure, NE, in the 43 GHz image, and 
detect a 0.75 mas jet feature in the 86 GHz image. The spectrum of the radio core is flat-to-steep 
and consistent with previous results, while the jet spectrum is flat-to-invert and inconsistent with 
one at lower frequencies. 
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Figure 2: Left : GMVA image of Mrk 501 at 86 GHz with model fit components with the restoring beam of 
205/i as x 37 uas in PA —5.11°. The peak brightness is 98 mJv beam 1 , and the contour levels are drawn at 
(-1, 1, 1.141, 2,aA?) x5.8 mJy beam 1 . The I a noise level is 1.9 mJy beam . Right : Spectral index map 
between the GMVA 86 GHz image and the third epoch VLBA 43 GHz image with total intensity contours 
of restored GMVA image. 

This observation cannot confirm the spectral index of the jet feature in Mrk 501 due to the 
limited (u, v)-coverage provided by the array used. To obtain more precise GMVA images, a better 
(u, v)-coverage, achievable by including new stations (e.g., the Korean VLBI Network) in the array, 
is indeed recommended. 
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